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We consider the paramagnetic phase of the random transverse-field Ising spin chain and study 
the dynamical properties by numerical methods and scaling considerations. We extend our previous 
work [Phys. Rev. B 57, 11404 (1998)] to new quantities, such as the non-linear susceptibility, higher 
excitations and the energy-density autocorrelation function. We show that in the Griffiths phase all 
the above quantities exhibit power-law singularities and the corresponding critical exponents, which 
vary with the distance from the critical point, can be related to the dynamical exponent z, the latter 
being the positive root of [(J//i)^''^] = 1. Particularly, whereas the average spin autocorrelation 
function in imaginary time decays as [G]av(T) ~ t"^''^, the average energy-density autocorrelations 
decay with another exponent as [G*]av(r) ~ t^^~^^^. 
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I. INTRODUCTION 



Quantum phase transitions occur at zero temperature 
by varying a parameter of the Hamiltonian, e.g. the 
strength of a transverse field. Quenched, i.e. time inde- 
pendent disorder has generally a profound effect on the 
properties of the quantum system not only at the critical 
point, but also in a whole region, which extends in both 
sides of the critical point. In this so called Griffiths phase 
the dynamical properties of the random quantum systems 
are exceptional: for example the (imaginary) time depen- 
dent average spin-spin correlations decay algebraically |^] 



(1.1) 



where the dynamical exponent z{S) is a continuous func- 
tion of the quantum control parameter, S. From here on 
we use [. . .]av to denote averaging over quenched disor- 
der. The physical origin of this type of singular behavior, 
as was pointed out by Grifhths Q for classical systems, 
is the existence of clusters in the random system, which 
are more strongly coupled than the average. The spins of 
such clusters, being locally in the "ordered phase", be- 
have coherently as a giant spin and the corresponding 
relaxation time is very large. Thus in an infinite system 
there is no (finite) time-scale and, as a con sequ ence, the 
autocorrelations decay algebraically, as in ( |l.l| ) j|] . 

Several physical quantities, which involve an integral 
of the autocorrelation function, (e.g. the static sus- 
ceptibility,) are singular not only at the critical point 
but also in a finite region of the paramagnetic phase. 
This phenomenon has been first noticed by McCoy in 
a two-dimensional classical model with correlated disor- 
der, (equivalent to a one-dimensional random quantum 
model) therefore we call the Grifhths singularities in 
quantum systems Griffiths-McCoy singularities. 



Many of the theoretical studies on random quantum 
systems are related to random quantum ferromagnets |^ 
and quantum spin glasses , which have also experimen- 
tal realizations In higher (d = 2 and rf = 3) dimen- 
sions one generally studies the distribution of the (linear 
and non-linear) susceptibilities, the asymptotic behavior 
of those can be related to the dynamical exponent, z{S), 
by scaling considerations. According to numerical stud- 
ies - in agreement with these phenonicnological theories - 
z{S) is found continuous function of the quantum control 
parameter S, which appear to have a finite limiting value 
at the critical point, S ~ 0, of spin glasses whereas it 
is diverging for random ferromagnets |^ . 

Many features of Grifhths-McCoy singularities can al- 
ready be seen in one-dimensional systems, where many 
exact and conjectured results exist. In this paper we 
consider the prototype of random quantum systems the 
random transverse-field Ising model (RTIM) in one di- 
mension, defined by the Hamiltonian: 



I 



I 



(1.2) 



Here the af , erf are Pauli matrices at site I and the J/ ex- 
change couplings and the hi transverse-fields are random 
variables with distributions 7r(J) and p{h), respectively. 
Note that in one dimension all the couplings and fields 
can be taken positive through a gauge transformation. 
The model in (1.2) is in the ferromagnetic (paramagnetic) 
phase if the couplings in average are stronger (weaker) 
than the transverse fields. As a convenient quantum 
control-parameter one can define: 

[ln/i]av-[lnJ] av 



<5 = 



(1.3) 



var[lnh] + var[ln J] ' 

and at the critical point 6 = 0. 

The Hamiltonian in eq( |l.2D is closely related to the 
transfer matrix of a classical two-dimensional layered 
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Ising model, which was first introduced and partially 
solved by McCoy and Wu js). Later the critical prop- 
erties of the quantum model was studied by Shankar and 
Murthy ||], and in great detail by Fisher [0. Through 
a renormalization group (RG) transformation Fisher has 
obtained many new results on static quantities and equal 
time correlations, which are claimed to be exact for large 
scales, i.e. in the vicinity of the critical point. Many 
of Fisher's results have been checked numerically pl| ] 
and in addition new results have been obtained about 
critical density profiles ||l^, time-dependent critical cor- 
relations [l^ ] and various probability distributions and 
scaling functions pl| , p^ . Later, using simple expressions 
about the surface magnetization and the energy gap sev- 
eral exact results have been derived by making use of a 
mathematical analogy with surviving random walks [ p^ , 
see also |ll5| . 

In the Griffiths phase, where the RG results are re- 
stricted to the immediate vicinity of the critical point, i.e. 
as (5 — > 0, numerical investigations both on temperature 
dependent [ p^ (specific heat, susceptibility) and dynam- 
ical quantities (spin-spin autocorrelations, distribution of 
the energy gap and susceptibility) ||l^,|llj have lead to the 
conclusion, that the behavior of all these quantities is a 
consequence of Griffiths-McCoy singularities and can be 
characterized by a single varying exponent z{S) in (LI). 
Very recently an analytical expression for z{6) has been 
derived by using an exact mapping [|l^ between the 
Hamiltonian in (1.2) and the Fokker-Planck operator of 
a random walk in a random environment. The dynami- 
cal exponent, which is given by the positive root of the 
equation: 



1/2 



(1.4) 



generally depends both on S and on the distributions 
7r(J) and p{J). However it becomes universal, i.e. distri- 
bution independent, in the vicinity of the critical point 
when z(5) w 1/(2(5), \6\ ^ 1, in accordance with the RG 
results Id. The numerical results obtained about dif- 
ferent singular quantities in the Griffiths phase are all in 
agreement with the analytical formula in ( |l.4| ) and the 
observed small deviations are attributed to finite-size cor- 
rections ll,!!. 

The singular quantities studied so far in the Griffiths 
phase are all related to the scaling properties of the low- 
est energy gap, which explains the observation why a 
single varying exponent is sufficient to characterize the 
singularities of the different quantities. There are, how- 
ever, other observables, which are expected to be singular 
too, but not connected directly to the first gap. For ex- 
ample one could consider the distribution of the second 
(or some higher) gap. By similar reasons as for the first 
gap these higher excitations are also expected to vanish 
in the thermodynamic limit and the corresponding prob- 
ability distributions are described by new exponents for 
small values of the gaps. As another example we consider 



the connected transverse spin autocorrelation function 
Gf(r) = ((gf ( 0)gf (r))). In the two-dimensional classi- 
cal version of ( |l.2[ ), the McCoy- Wu model, this function 
corresponds to the energy-density correlation function in 
the direction, where the disorder is correlated. There- 
fore we adopt in the following this terminology and call 
G^(r) the energy-density autocorrelation function. Since 
the inverse time scale for these correlations is, as we shall 
see, determined by the second gap, one expects that also 
[G'^]av(''') has an algebraic decay: 



(1.5) 



with an exponent rj^. Finally one should mention the 
non-linear susceptibility whose distribution is expected 
to be described by a new varying exponent. 

In this paper we extend previous numerical work and 
study the scaling behavior of the above mentioned singu- 
lar quantities in the Grifhths phase. We present a phe- 
nomenological scaling theory and we confront its predic- 
tions by results of numerical calculations, based on the 
free-fermion representation of the Hamiltonian in ( |1.2| ). 
We show that the physical quantities we studied are char- 
acterized by power-law singularities with varying critical 
exponents, the value of those are connected to the dy- 
namical exponent through scaling relations. 

Throughout the paper we use two types of random 
distributions. In the symmetric binary distribution the 
couplings could take two values A > 1 and 1/A with the 
same probability, while the transverse-field is constant: 



n{J)^^{SiJ-X)+SiJ-X-')) , 
p{h) - 6{h - ho) 



(1.6) 



At the critical point ho — I, whereas in the Griffiths 



phase, 1 < ho < X, the dynamical exponent from (1.4) is 
determined by the equation: 



ho^^ = cosh 



In A 



(1.7) 



In the uniform distribution both the couplings and the 
fields have rectangular distributions: 



pih) 



h. 



1, for < J < 1 
0, otherwise 

\ for < h < h 
0, otherwise 







(1.8) 



The critical point is also at /iq = 1, whereas the dynam- 
ical exponent is given by the solution of the equation: 

zln{l - z-^) = -Inho , (1.9) 

where the Griffiths phase now extends to 1 < /lo < cxd. 

The structure of the paper is the following. In Sec- 
tion 2. we present the free fermion description of various 
dynamical quantities. Phenomenological and scaling con- 
siderations are given in Section 3, the numerical results 
are presented in Section 4. Finally, we close the paper 
with a Discussion. 
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II. FREE FERMION DESCRIPTION OF 
DYNAMICAL QUANTITIES 

We con sider the random transverse-field Ising model 
in eq(1.2) on a finite chain of length L with free bound- 
ary conditions. The Hamiltonian in eq(1.2) is mapped 
through a Jordan- Wigner transformation and a following 
canonical transformation |19f| into a free fermion model: 



9=1 



(2.1) 



in terms of the rf^ (rjq) fermion creation (annihilation) op- 
erators. The energy of modes, eg, is obtained through the 
solution of an eigenvalue problem, which necessitates the 
diagonalization of a 2L x 2L tridiagonal matrix with non- 
vanishing matrix-elements T2i-i^2i = '72i,2i-i = hi, i — 
1, 2, . . . , L and T2i^2i+i = T2i+iai = « = 1, 2, . . . , L— 1, 
and denote the components of the eigenvectors Vq as 

Vq{2i - 1) = and Vq{2i) = i = 1,2,...,L, 

i.e. 





( 


hi 


\ 




hi 





■h 







Ji 


h2 


T = 






h2 








Jl-i 








Jl-1 hL 



/ 



-*9(2) 



Va 



\ 



(2.2) 



-%{L) 



We consider only the > part of the spectrum |2C|] . 

The local susceptibility xi at site I is defined through 
the local magnetization m; as: 



XI 



hm — — 

Hi^o 5 Hi 



(2.3) 



where Hi is the strength of the local longitudinal field, 
which enters the Hamiltonian (l.S ) via an additional term 
Hiaf . xi can be expressed as: 



XI 



Ej — En 



(2.4) 



where |0) and \i) de not e the ground state and the ith ex- 
cited state of -ff in (1.2) with energies Eq and Ei, respec- 
tively. For boundary spins one has the simple expression: 



(2.5) 



Similarly, the local non-linear susceptibility is defined by: 

S^mi 



xf = lim , , 
H,^o 6Hf 



(2.6) 



and can be expressed as: 



It should be noted that it is not the first sum on the r.h.s. of ( ^.7[ ) that gives the leading contribution, since at leas t 
one of the 3 energy differences most involve a higher excitation {{i\(j'f\i) = for i = j). For surface spins, I = 1, (2.7) 
simplifies: 



(2.8) 



Next we consider the energy- density correlation function = ^ KOjcrf |0)p exp[— r(i?i — Eq)] . (2.9) 

at site Z, G^, defined by: i>o 

Gf(T) = (0|CTf (T)crf (0)10) - (0|(Tf (T)|0)(0|crf (0)|0) In the free-fermion representation it is given by: 
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i5>7 

(2.10) 

which can be expressed for surface spins as: 



<5>7 



£5 



hi 



-05(1)07(0 



exp[-r(e5 + e^)] 



(2.11) 



The spin-spin a utoc orrelation function, Gi, which is de- 
fined as Gf in (2.9) by replacing af by erf, is generally 
complicated and can be expressed in the form of Pfaffi- 
ans ||2l|,|l^ . Exception is the autocorrelation function for 
surface spins, which is simply given by: 



Gi(T)-5]|$,(l)pexp(-re,). 



(2.12) 



III. PHENOMENOLOGICAL AND SCALING 
CONSIDERATIONS 

As described in the Introduction the Grifhths-McCoy 
singularities in the paramagnetic phase are connected 
to the presence of strongly coupled clusters, which are 
locally in the "ordered phase" and therefore the cor- 
responding excitation energy is very small. For the 
RTIM the origin of these clusters can be explained either 
through the analysis of the RG fixed-point distribution 
pO| , which works only in the vicinity of the critical point, 
or by using simple explicit expressions for the excitation 
energy p^, |l4| and estimate those through random walk 
Here we use a simple phenomenological 
23| ], whose results are in agreement with 
the above microscopic methods. 

Consider the quantity, Pl{N), which measures the 
probability that in a chain of L sites there is a cluster 
of <C L strongly coupled spins. Since N consecutive 
strong bonds can be found with exponentially small prob- 
abihty ~ exp(— AiV), whereas the cluster could be placed 
a.t ^ L different sites we have 



arguments ]14 
approach [ll[' " 




Pl(A^) - iexp(-yliV) 



(3.1) 



The excitation energy of this sample corresponds to the 
energy needed to flip all spins in the cluster, which is 
exponentially small in N: 



ei ~ exp{-BN) 



(3.2) 



Combining (3.1) with ( |3.2| ) we have for the probability 
distribution of the first gap 



Pi(lnei) ^Lei 



l/z 



(3.3) 



for ei ^ and l/z = A/B. Here, from the scaling com- 



bination in eq(3.3): L ^ ^ r^/^, we can identify z 

as the dynamical exponent. 

Next, we consider the second gap, €2, which is con- 
nected to the existence of a second strongly connected 
cluster ol N' < N spins, and its value corresponds to the 
energy needed to flip all the spins in the second cluster 
simultaneously, consequently 



62 ~ exp(-SiV') 



(3.4) 



The probability with which a cluster of size N' occurs, 
provided another cluster of size N > N' exists, is given 

by PiiN')^^ iexp(-A7V') ELjv' PlW. For N' « L 
(or in the infinite system size limit L 00) this can be 
estimated as: 



P^(iV') - L'^ exp[-2AN'] 



Thus from (pA) and (3.5) we have 



P^lnes) 



(3.5) 



(3.6) 



(3.7) 



with l/z' = 2A/B, thus 

z' = z/2 



Note that the scaling combination in the r.h.s. of ( |3.6| ) 
is dimensionless, as it should be. Repeating the above 
argument for the third, or generally the nth gap the 
corresponding distribution is described by an exponent 
^(n) _ ^y'j^^ however the finite size corrections for these 
gaps are expected to increase rapidly with n. 

The scaling behavior of the probability distribution of 
the susceptibilities can be obtained by noticing that both 
for xi E^nd x"' the leading size dependen ce i s conn ecte d 
with energy gaps in the numerator of (2^) and (^), 
respectively. Then for the asymptotic behavior of the 
distribution of the local susceptibility we have: 



1 



In [P(lnx/)] = InXi + const , 



(3.8) 



like to the inverse gap. For the no n-li near susceptibility 
the second term in the r.h.s. of Eq.(2/7) gives the singular 
contribution, so that 



ln[P(lnxr')] 



1 



„i iT^xf + const , 



with 



= 3z 



(3.9) 



(3.10) 



since the asymptotic distribution is the same as that of 
the third power of the inverse gap. We note that the 
relation in (3.1C) corresponds to the phenomenological 
result in ^. 

The scaling behavior of the average spin autocorrela- 
tion function is given by: 
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[G;]av(T) 



yPL(ei)|M,pexp(-Tei)dei , (3.11) 



where the factor with the matrix-element is |M;p ~ 1/L, 
since the probability that a low energy cluster is local- 
ized at a given site, is inverse ly p roportional to the 
length of the chain. Then using ( |3.3| ) one arrives to the 
result in (LI), thus establishing the relation between the 
decay exponent of the spin autocorrelation function and 
the dynamical exponent. 

For ene rgy-d ensity autocorrelations, according to 
.10) and (2.11) the characteristic energy scale is £2 and 
the asymptotic behavior of the average energy-density 
autocorrelation function is given by: 



(1 



J Pae2)|Mf|2exp(-re2)de2 ■ (3.12) 



Now we take the example of the surface autocorrelation 
function in (2.11) to show that the factor with the matrix- 
element, \ AIf is proportional to e^- The remaining 
factor in (2.11) with the first components of the eigen- 
vectors is expected to scale as 1 /L due to similar reasons 
as for the spin autocorrelations, thus {M^l"^ ~ L~^e2 and 



together with (|3J) one has (e2)|Mf ^ - Le 



l/z' + l 



Be- 
fore evaluating the integr al in (3.12) we note that for a 
fixed L the expression in (3.12) stays valid up to t ~ L^. 
Therefore to obtain the L independent asymptotic be- 
havior in T we should instead vary L, so that according 

to (3.6) take L ^ €2 ^^^^^ and in this way we stay within 
the border of validity of ( 3. IS ) for any r. With this mod- 
ification we arrive to the result in Eq.(1.5) where the 
decay exponent, rje, is related to the dynamical exponent 
as 



z 



(3.13) 



where the relation in (3.7) is used. We expect that the 
factor, \Mf\^, has the same type o f sca ling behavior for 



any position I, thus the relation in ( ^.13 ) stays valid both 
for bulk a nd su rface spins. We note that the reason- 
ing above (3.13) applies also for the spin autocorrelation 
function, in which case in ( ^.ll| ), however there is no ex- 
plicit L dependence. 

In this way we have established a phenomenological, 
scaling theory which makes a connection between the un- 
conventional exponents in (^), ( p. 10 ) (3.13) and the dy- 
namical exponent. In the next Section we confront these 
relations with numerical results. 



We start to present results about the distribution of 
the first and second gaps. As illustrated in Fig. 1, both 
for the uniform and the binary distributions, the asymp- 
totic sca ling relations for the distribution of the first two 
gaps in (3.3) and (3.6) are satisfied. From the asymptotic 
slopes of the distributions in Fig. 1 we have estimated 
the 1/z and 1/z' exponents for the two largest finite sys- 
tems, L = 64 and L = 128, which are presented in Fig. 
2 for different points of the Griffiths phase for the uni- 
form distribution. As seen in the Figure the z exponent 
calculated from the first gap agrees very well with the 
analytical results in (l.£). For the z' exponent, as calcu- 
lated from the distribution of the second gap the scaling 
result in ( |3.7[ ) is also well satisfied, although the finite- 
size corrections are stronger than for the first gap. For 
the third gap, due to the even stronger finite size effects, 
we have not made a detailed investigation. Extrapolated 
results at /lo = 2 are found to follow the scaling result 
^(3) = z/3. 

Next, we study distribution of the linear and non-linear 
local susceptibilities at the surface spin. As demonstrated 
in Fig. 3 both type of distributions sat isfy the respective 
asymptotic relations in ( |3.8| ) and (3.9), from which the 
critical exponents z and z"' are calculated. The esti- 
mates are shown in Fig. 4 at different points of the Grif- 
fiths phase. As seen in the Figure the numerical results 
for the dynamical exponent, z, are again in very good 
agreement with the analytical results in (1.9) and also 
the exponent of the non- linear su scept ibility, z"', follows 
fairly well the scaling relation in ( 3.10|) . 



Finally, we calculate the average energy-density auto- 
correlation function. As seen in Fig. 5 [G"^]av(''') displays 
a linear region in a log-log plot, the size of which is in- 
creasing with L, but its slope, which is just the decay 
exponent, jye, has only a weak L dependence. The slope 
of the curve and thus the corresponding decay exponent 
r\e has a variation with the parameter /ig, as illustrated 
in Fig. 6. The estimated Tye exponents at the critical 
point, /iQ = 1, and in the Griffiths phase are presented in 
Fig. 7. As seen in this Figure the variation of r/e is well 
described by the form 7ye(<5) — 7?e(0) + l/z((5). This func- 



tional form corresponds to the scaling result in Eq.(3.13) 
if the critical point correlations decay with 



r]M = 2 . 



(4.1) 



IV. NUMERICAL RESULTS 

In the numerical calculations we have considered RTIM 
chains with up to L = 128 sites and the average is per- 
formed over several 10000 realizations, typically we con- 
sidered 50000 samples. For some cases, where the fi- 
nite size corrections were strong, we also made runs with 
L = 256, but with somewhat less realizations. 



The numerical calculations with L = 128 give a slightly 
higher value rye(O) ~ 2.2 |Q. However, the finite-size 
estimates show a slowly decreasing r]e{0) with increasing 
system size. Repeating the calculation with L = 256 we 
got r?e(0) « 2.1. Thus we can conclude that the scal- 



ing relation in (3.13) is probably valid and then Eq.(4.1) 
is the exact value of the decay exponent of the average 
critical energy-density autocorrelations [e3. 
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V. DISCUSSION 



In this paper we have considered the random 
transverse-field Ising spin chain and studied different 
consequences of the Griffiths-McCoy singularities in the 
paramagnetic phase. Our main conclusion is that all sin- 
gular quantities can be characterized by power-law sin- 
gularities and the corresponding varying critical expo- 
nents can be related to the z{S) dynamical exponent and, 
for energy-density autocorrelations, to the ??e(0) critical 
poi nt e xponent. Since the exact value of z(5) is known 
in ( |l.4|) and we expect that also the relation in (O) is 
valid, thus wc have a complete, analytical description of 
the Griffiths phase of the RTIM in one dimension. 

One interesting feature of our results concerns the dis- 
tribution of the higher excitations and the value of the 
corresponding exponent z^"-* — z/n. Since the decay of 
dynamical correlations of general, more complex opera- 
tors are related to l/z^"^ = n/z, we obtain a hierarchy 
of decay exponents which could be simply expressed by 
those of a few primary operators. This feature is reminis- 
cent to the tower-like structure of anomalous dimensions 
in two-dimensional conformal models [^5|| . Our knowl- 
edge about the higher excitations can also be used to 
estimate the correction to scaling contributions. 

Much of the reasoning of our phenomenological scaling 
considerations in Section 3 stay valid for other random 
quantum systems. Especially the scaling behavior of the 
higher gaps and the corresponding relation in Eq.(3.7) 
should be valid even for higher dimensions and the same 
is true for the distribution of the no n-lin ear susceptibility 
and the corresponding relation in ( [3.10[ ). 

In one dimensions the universality class of the RTIM 
involves several random systems, among others the ran- 
dom quantum Potts chain |2^. For these models one 
docs not expect a universality of the z{d) exponent in 
the Grif fith s phase, although scaling relations, like the 
one in (3.7) are very probably valid. It would be in- 
teresting to perform a numerical study on the random 
quantum Potts model to check the existing conjectures. 

Another possible field where the present results could 
be applied is the problem of anomalous diffusion in a 
random environmen t [ r7in^ ,|27|] . Making use of the ex- 
act correspondence |p7,|18| between the Hamilton opera- 
tor in (|1.2|) and the Fokker-Planck operator of the one- 
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dimensional random walk we can use the relation in (3.7) 
to describe the distribution of the eigenvalues of the cor- 
responding Fokker-Planck operator. One can also define 
analogous quantity to the energy-density autocorrelation 
function in ( 2.10| ) by considering connected persistence 
correlations, whose asymptotic decay i s rel ated to the 
distribution of the second gap, like in ( 3.12| ). Research 
in this field is in progress. 
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FIG. 2. The estimates for 1/z and 1/z' as a function of 
ho for the uniform distribution. They have been obtained 
from our analysis of the probability distribution of Inei and 
In £2, respectively, for two system sizes (as exemplified in fig. 
1). The full line for 1/z corresponds to the analytical result 
(1.4), the broken line corresponds to 2/z, which we predict to 



be identical with 1/z' . 
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FIG. 3. Probability distribution of the linear and 
non-linear susceptibility, In^i and Inxi', respectively, for the 
uniform distribution at ho = 3. The straight lines are least 
square fits to the data for the largest system size, their slopes 
correspond to l/z{ho) and l/z"'(/io), respectively. They fol- 
low the predicted relation «"'(/io) = 3z{ho). 
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FIG. 1. Probability distribution of Inei and In £2 for the 
uniform distribution nt ho — 2 (top) and the binary distri- 
bution (A = 4) at ho — 2.5 (bottom). The straight lines are 
least square fits to the data for the largest system size, their 
slopes correspond to l/z{ho) and 1/ z' {ho), respectively. They 
follow the predicted relation z' {ho) = z{ho)/2. 
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FIG. 4. The estimates for 1/z and l/z"' as a function 
of ho for the uniform distribution. They have been obtained 
from our analysis of the probability distribution of In xi and 
Inxi', respectively, for two system sizes (as exemplified in fig. 
3). The full line for 1/z corresponds to the analytical result 
(1.4), the broken line corresponds to l/Sz, which should be 



identical with 
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FIG. 5. The bulk energy-energy autocorrelation function 
[GL/2lav(T) for the binary distribution (A = 4) at ho — 1.5 
for different system sizes as a function of Inr. The slope of 
the straight line identifies the exponent rj^{ho) describing the 
asymptotic decay of [G2/2]av(T) in the infinite system size 
limit L —f cx>. 
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FIG. 6. The bulk energy-energy autocorrelation function 
[G2/2lav(T) for the binary distribution (A = 4) at different 
values for ho for L = 128 as a function of Inr. One observes 
the variation of the exponent rje.(ho) (identical to the slope of 
the straight line fits) with increasing ho- 
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FIG. 7. The exponent rjeiho) for the binary distribu- 
tion (A = 4) as obtained from the analysis of the asymp- 
totic decay of the bulk energy-energy autocorrelation func- 
tion [G2/2]av(T) a la fig. 6. The full line is the analytical 
predictio n r)^ (ho) — 2 + l/z{ho) with z{ho) given by the exact 
formula (1.4) for the binary distribution with A = 4. 
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